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The recently assembled laser-beam detectors of gravitational waves are approaching the planned 
level of sensitivity. In the coming 1-2 years, we may be observing the rare but powerful events of 
inspiral and merger of binary stellar-mass black holes. More likely, we will have to wait for a few 
years longer, until the advanced detectors become operational. Their sensitivity will be sufficient 
to meet the most cautious evaluations of the strength and event rates of astrophysical sources of 
gravitational waves. The experimental and theoretical work related to the space-based laser-beam 
detectors is also actively pursued. The current gravitational wave research is broad and interesting. 
Experimental innovations, source modelling, methods of data analysis, theoretical issues of principle 
are being studied and developed at the same time. The race for direct detection of relatively 
high-frequency waves is accompanied by vigorous efforts to discover the very low-frequency relic 
gravitational waves through the measurements of the cosmic microwave background radiation. In 
this update, we will touch upon each of these directions of research, including the recent data from 
^ , the Wilkinson Microwave Anisotropy Probe (WMAP). 
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I. INTRODUCTION 



CN| ■ The concept of gravitational radiation has been with us for quite a long time. Thinking of the relativistic gravita- 
tional field in parallel with the familiar case of the electromagnetic field, it was natural to expect that there should 
exist waves of the gravitational field similar to the waves of the electromagnetic field. As Einstein [1] put it in 1913: 
"The conviction had to come that Newton's law of gravitation is as incapable of desribing all gravitational phenomena 
as Coulomb's laws of electrostatics and magnetostatics are of electromagnetic phenomena". The decades of hard work 
have followed. In the beginning, the research was purely theoretical. Through doubts and controversies, the concep- 
ts") , tual and mathematical issues of gravitational radiation have been clarified. Then, in the 60-s, the experimental work 
has started with the pioneering effort of J. Weber [2]. In the 90-s, gravitational waves have been observed indirectly 
""^ | via the measurement of secular changes in the orbital parameters of the binary system of neutron stars that includes 
Q H . the pulsar PSR 1913+16 [3]. These days, gravitational waves are routinely taken into account in the theoretical and 
observational studies ranging from orbital evolution of close pairs of compact stars to the early Universe cosmology. 
The experimental progress has also been very impressive. 

We are now at a special and decisive p[oint. First, the relevance and importance of gravitational- wave research is 
fully recognised by the communities of physicists and astronomers. The worldwide network of scientific collaborations 
has been established, reflecting the necessity of coincident observations at various instruments and the need for 
joint analysis of the data [4], [5]. Second, the scientific runs have begun at the recently assembled sensitive laser- 
interferometric observatories - American LIGO [6] and British-German GEO600 [7]. The French-Italian VIRGO [8] 
will be operational soon. Meanwhile, the Japanese TAMA300 [9], along with the international network of bar detectors 
[10], continue to collect data at their level of sensitivity. The design sensitivity of three LIGO interferometers plus 
VIRGO and GEO600 meets the realistic astrophysical predictions, so the direct detection of powerful (even if rare) 
sources, such as coalescing stellar-mass black holes, becomes likely. Third, the European Space Agency (ESA) and 
American NASA have agreed to share the costs of the space-based Laser Interferometer Space Antena (LISA) [11]. 
LISA is planned to be launched around the year 2011, preceeded by a technology demonstration mission. The plans 
for the advanced ground-based detectors, such as LIGO II, the Japanese cryogenic laser interferometer LCGT, and, 
possibly, the European EURO, are also maturing very quickly. These instruments of the next generation will detect 
a host of well anticipated sources, such as compact binary stars, but really fundamental discoveries are also expected. 
Fourth, there exists a strong competition on the side of purely astronomical means of indirect detection of gravitational 
waves. In fact, the anisotropy and polarisation measurements of the cosmic microwave background radiation (CMB) 
are likely to bring us decisive information on the fundamentally important relic gravitational waves much earlier than 
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it will be done by direct methods. By all counts, gravitational-wave research is now one of the most exciting and 
promising areas of physical science. 

The weakness of gravity as a physical interaction is the strength of gravitational waves as a tool of scientific 
research. It is difficult to detect gravitational waves because they carry their energy practically without scattering 
and absorption. But this is exactly because of this difficulty that we have any chance to learn about what was 
happening at the beginning of the Universe, or in the depth of an exploding supernova, or in the vicinity of, what 
we think are, merging black holes. It is also likely that the differing properties of gravitational waves will be a 
discriminating signature of different fundamental physical theories and modified gravities. 

In this review, we will start, in Sec. 2, with an elementary theory of gravitational waves. Then, in Sec. 3, we 
will discuss the current status of gravitational-wave experiments. In Sec. 4 we will focus on astrophysical sources of 
gravitational waves and the new physics that will be learned from their observation. Sec. 5 is devoted to gravita- 
tional waves and cosmology. In particular, some recent results from the WMAP will be discussed in the context of 
gravitational-wave research, including the "implications for inflation" . Finally, we briefly summarise this update in 
Sec. 6. The background material on the gravitational-wave research can be found in [12], [13], [14], and some of the 
previous reviews in [15], [16], [17], [18], [19]. 

II. ELEMENTARY THEORY OF GRAVITATIONAL WAVES 

There are some similarities in the mathematical description of the gravitational field and the electromagnetic field. 
The electromagnetic field can be described by the components of the 4- vector potential A^. The quantities are 
functions of time and spatial coordinates x a , x a = (ct, x, y, z), and they obey the wave-like dynamical equations - the 
Maxwell equations. Other quantities, such as electric and magnetic components E, H, and the energy-momentum 
tensor of the electromagnetic field, are calculable from A^ . The quantities A 11 allow the gauge freedom. This 
means that some seemingly different solutions of the field equations are in fact equivalent solutions in the sense of 
their physical manifestations. Analogously, the gravitational field can be described by 10 components of the 4x4 
symmetric tensor as functions of x a . The quantities h tiV obey the nonlinear wave- like dynamical equations - the 
Einstein equations. The energy-momentum tensor of the gravitational field #*" is calculable from the tensor h^ u . The 
quantities W allow the gauge freedom, so that some seemingly different solutions of the gravitational field equations 
describe in fact the physically equivalent configurations. 

In contrast to the electromagnetic field, the gravitational field W is a nonlinear field. This means that a sum of 
two solutions of the field equations is not a new solution, and the gravitational field, along with matter fields, is a 
source for itself. However, in a number of situations one can neglect the nonlinearity of the gravitational field. In this 
approximation we come to the notion of weak gravitational fields and linearised gravitational waves. 

Many properties of linearised gravitational waves resemble those of electromagnetic waves. Gravitational waves 
propagate with the velocity of light c and have two independent transvers polarization states. In its action on 
free masses, a gravitational wave (g.w.) exhibits some analogs of the electric and magnetic contributions of an 
electromagnetic wave (em.w.) acting on free electric charges. Gravitational wave field is dimensionless and its strength 
can be characterized by a dimensionless amplitude h. The amplitude h decreases in the course of propagation from 
a localized source in inverse proportion to the traveled distance: h oc 1/r. Gravitational waves carry away from a 
radiating system its energy, angular momentum and linear momentum. 

The linearised g.w. satisfy the wave equation 

h ^« a + ^ h ^ a (i _ h ^ a _ h ^ a = o, (1) 
where the ordinary derivative is denoted by a comma and jf" is the metric tensor of the Minkowski space-time: 

dcr 2 = rj^dx^dx" = c 2 dt 2 - dx 2 - dy 2 - dz 2 . (2) 
The first term in Eq. (1) is the familiar d'Alembert (wave) operator. A plane-wave solution to Eq. (1) is given by 

hf™ = a' M/ e ikaXa , (3) 

where k a k a — 0, reflecting the fact that a g.w. propagates with the velocity of light. Because of this condition, 
the field equations (1) require the 10 components of the constant matrix a^ u to satisfy 4 constraints: a^ v k v = 0. 
Then, the quantities a 00 and a 0z can be expressed in terms of 6 components of the matrix a 1 - 7 '. This matrix itself 
includes the part a n satisfying the further 4 constraints: a lJ kj = 0, a lJ riij = 0. These remaining 2 degrees of freedom 
(sometimes called the TT-componcnts) fully determine the observational manifestations of the plane wave and its 
energy-momentum characteristics. Indeed, it is easy to show that the gravitational energy-momentum tensor 
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depends only on the TT-components of the field: 
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where we have dropped (as we normally do in the case of electromagnetic waves) the purely oscillatory terms. 

How does a gravitational wave affect the free (i. e. not subject to any other forces or constraints) masses ? In 
the case of an em.w., we could have noticed the displacement of a charged particle with respect to a neutral particle 
placed initially in the same point. In the case of a gravitational wave, there is no particles neutral to the gravitational 
interaction, so we need to study the relative displacement of particles separated initially. This brings us to the analysis 
of tidal effects of a gravitational wave, similar to the tidal effects of a Newtonian gravitational field. Let one of free 
masses define the origin of our coordinate system, xfa = 0, whereas the second mass is placed initially at x\ 2 \ = l % 

and its (small) displacement caused by the gravitational wave is denoted by £*. Then the equations of motion of the 
second mass are: 



4^ = ~u,Vj fe e iwt . 
dt 2 2 k 



(6) 



These equations depend only on the TT-components of the wave. 

For a wave propagating, say, in z-direction it is convenient to write the TT-components explicitely: 



"vv 



h + sin(wi — kz + -0), h xy — h x cos(wi — kz + tp) 



(7) 



Then, the relevant solution to Eq. (6) reads 



x = l\ — 2^1^+ si 11 ^ + "0) — 2^ 2 ^ x cos ( LlJ ^ + VOj V = h 



—lih x cos(ujt + ijj) + —l 2 h + sm(ujt + ip), z — l 3 . (8) 




FIG. 1. Motion of free particles in the field of a linearly-polarised gravitational wave 



The masses that lie (on average) on the ring (^i) 2 + (I2) 2 = I 2 , h = 0, oscillate around their ring positions. In Fig. 
1 we show a quarter of the cycle caused by a wave with the + polarization state (i.e. when h + ^= 0, h x = 0). The 
pattern of oscillations enforced by the x polarization state (i.e. when h x 7^ 0, h+ = 0) can be obtained from Fig. 1 
by its rotation by the angle 45° in the (x, y) plane. The general motion of a particle is a linear superposition of these 
oscillations, as seen in Eq. (8). The oscillatory deformation of a sphere of masses, surrounding the central mass at 
the origin, is desribed by a liner combination of spherical harmonics Yj m (#, (f>) with I = 2 and m = ±2, and where the 
polar axis is taken along the z direction. 

Equations (8) suggest that the motion of free particles is confined strictly to the planes z — const. This conclusion 
comes about only because we have so far neglected the smaller terms, containing the products oih+, h x with the extra 
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small factor l/X, where / is the separation between masses and A = 2ttc/uj is the gravitational wavelength. When these 
terms are taken into account, the perturbed positions x, y and, most importantly, z receive further contributions. The 
gravitational wave drives the masses not only in the plane of the wave-front, but also, to a smaller extent, back and 
forth in the propagation direction [16]. This extra component of motion is similar to the one caused by the magnetic 
field and the Lorentz force of the em.w. acting on a charged particle. In Fig. 2 we show typical displacements of 
the masses forming (on average) the plane z = 0, when the "magnetic" contribution to their motion is also taken 
into account. It is seen from this figure that the entire plane of particles is being bent in an oscillatory fashion. This 
"magnetic" component of motion will play a certain role in the analysis of observations with laser interferometers [20] . 
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FIG. 2. Motion of free particles with the "magnetic" component taken into account 

The amplitudes h+, h x are determined by the source of the gravitational waves. To find the amplitudes, we should 
replace the zero in the right hand side of Eq. (1) by the source term (167rG/c 4 )T A " y , where T^ v is the energy- 
momentum tensor of the source, and seek the retarded solutions to the wave equation. Assuming that the distance 
i?o to the source is much larger than the wavelength, one can write 

i£tf"<*.) -*"<*.)]. 

where the Fourier components T^ v (k a ) arc defined by 



h+(k a ) = — — [T 11 (k a )-f 22 (k a )}, h x (k a ) = ——T 12 (k a ), (9) 



T^(k a )= J (J T^(t r ,r )d 3 r Q je- lk «* a d 4 x 

and t r is retarded time. These formulas allow one to evaluate the typical amplitude h from a given source: 

1 GM fv\ 2 



h . 



Ro c 2 V c > 

where M is the total mass of the source and v is the characteristic non-spherical velocity of the matter bulk motion. 
For a bound system like a binary star, (v/c) 2 ~ GM/c 2 a, where a is the size of the system. An accurate calculation 
for a binary in a circular orbit, consisting of masses Mi, M 2 separated by the distance a, and after averaging over the 
orbital period and orientation of the orbital plane, gives 

, /n2» ,,2^1/2 /32\ 1/2 1 G 5 / 3 M!M 2 . ,. 2/3 

where the emitted g.w. frequency / (in Hz) is 

/ = i 

Note that the amplitude h depends on a particular combination of masses: 

MlM " -M*'\ where M = ,^ M f* (13) 



G{M 1+ M 2 )Y'\ (12) 



(Mi+M 2 )!/3 ' (Mi+M 2 )!/5' 

so that h oc (l/Ro)M 5 / 3 f 2 / 3 and M. is sometimes called a chirp mass. 
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III. CURRENT STATUS OF GRAVITATIONAL WAVE DETECTORS 



The possible methods of detecting gravitational waves, as well as the associated difficulties, can be seen from the 
discussion above. Fig.l is helpful for understanding the principles of mechanically coupled detectors (bar detectors) 
and electromagnetically coupled detectors (laser interferometers). A bar detector is essentially a mechanical oscillator 
consisting of two masses connected by a spring. As an illustration, one can think of two elastically connected masses 
lying, say, at the opposite ends of the x axis in Fig. 1. The size of bar detectors is normally small, a couple of 
meters or so. A bar detector is a relatively narrow-band instrument. It is mostly sensitive to g.w. frequencies in 
the vicinity of the main eigen-frequency of the bar. In the presently operating instruments the resonant frequency is 
around ~ 1kHz. There are some advantages in using spherically shaped mechanical detectors. The construction of 
such detectors is currently taking place (see, for example, [21]). While the operating bar detectors continue to collect 
useful information [22], [23], [24], we will concentrate on laser interferometers. 

The laser interferometer technique is based on free masses- mirrors whose relative distances are monitored by bounc- 
ing light. A ground-based interferometer is normally an L-shaped configuration. The light beamsplitter is in the corner 
of L, while the reflecting mirrors are at the ends of the configuration and near the beam splitter. The central mass 
in Fig.l is an illustration of the beamsplitter and corner mirrors, whereas the end mirrors are located, say, in the 
positive directions of x and y. Obviously, mirrors in real interferometers are not free masses, they are suspended like 
pendulums. But they behave essentially as free masses in their motion along the corresponding arm. The length of 
the arms of an interferometer can be large. For instance, it ranges from 600 meters in GEO600 to 4km in largest of 
LIGO interferometers. In contrast to bar detectors, laser interferometers are relatively broad-band instruments. They 
are sensitive to frequencies in the interval w (30 — 10 4 )Hz. Fig.l shows the simplest case of a single monochromatic 
wave with one polarisation state arriving from the orthogonal direction to the detector's plane, but the response of 
the detector to the general case of the incoming wave is also calculable. 

It is clear from Eq. (8) that the relative displacement of free masses is proportional to the incoming wave amplitude 
h: 51/1 w h. Which numerical values of h can be expected from astrophysical sources in the surrounding Universe ? 
Consider one of the most powerful and efficient emmitters - a pair of compact stars orbiting each other at a tight 
orbit. Specifically, consider two neutron stars with masses M\ — M 2 = 1.5M at the late stage of their inspiral. 
Orbiting each other at separation a = lOOfcm, they emit gravitational waves at frequency / = 200Hz. The emitted 
intensity of radiation is very high by astronomical standards: 3 x 10 52 erg/ sec. Let the distance to the source be 
i?o = 100-Mpc. We cannot take this distance much shorter than 100-Mpc, because the expected number of such events 
per year in a smaller volume of the surrounding Universe would be less than 1. Then, formula (11) says that the 
amplitude at Earth is h w 10~ 22 . This is an increadibly small number. It enters any conceivable method of detection 
of gravitational waves and explains why it is so difficult to observe them. In a 4km long interferometer, we need to 
beat all the noises and measure the mirror's displacements at the level of 4 x 10 _ 17 cm. 

The interferometer monitors the time-dependent difference of the relative distance variations in the two arms. Using 
the terminology of solid state physics, this quantity is called the (dimensionlcss) strain. Regardless of the presence 
or absense of the useful g.w. signal, there will always be some strain noise n(t) in the output of the detector. The 
mean-square value of the noise can be expressed as an integral over the noise power spectral density S n (f): 

/>oo 

n 2 (t) = 2 / S n {f)df. (14) 
Jo 

The square root of S n (f) is the noise amplitude y/ S n (f). This quantity has the dimensionality of -ffz~ 1//2 , and it is 
this quantity that is usually plotted on the sensitivity graphs. In Fig. 3 (taken from [6]) we show the recent status of 
the noise amplitudes in the 4km interferometer of the LIGO Livingston cite. One can see the great progress that has 
been made on the way of reaching the design goal of sensitivity, shown by the solid line. The lower frequency part of 
the sensitivity curve is dominated by seismic perturbations, the central part by thermal noise in mirror's suspensions, 
and the higher frequency part by the shot noise of the laser light. The advanced LIGO sensitivity curve will be lower 
than the solid line In Fig. 3 by a factor of 10 across all the frequencies. This upgrade of LIGO is planned to take place 
in about the year 2007. The advanced instruments will also allow optical configurations in which the sensitivity can 
be increased in a certain narrow frequency band at the expense of lowering the sensitivity outside the chosen band. 

To compare qualitatively the astrophysical signal S, represented by the dimensionless amplit ude h, with the de- 
tector's r.m.s. noise N, we should calculate, as formula (14) suggests, the product J S n (f)\/ Af , where Af is 
the appropriate bandwidth in the noise sp ectrum . Taking the initial LIGO design sensitivity 3 x 10~ 23 .ffz~ 1//2 at 
/ = 200Hz and Af = f, we would get y/S n (f)y/A~J = N = 5 x 10~ 22 . This N is a factor of 5 higher than the 
signal S = 10~ 22 from the coalescing neutron stars considered above. The ratio S/N —1/5 would suggest that, 
by a significant margin, the signal is not measurable. The reality, however, is somewhat better than this estimate. 
One should take into account the apriori knowledge of the expected waveform. If the waveform is known, one can 
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use the well developed technique of matched filtering. This method recovers the signal by, effectively, reducing the 
appropriate A/ and the relevant amount of detector's noise. For a quasi-periodic signal, this is achieved by using a 
long observation time T: 1/T ~ A/ <C /. In the case of a coalescing binary, the dominant g.w. frequency is increasing 
with time, but the binary still executes almost n = 200 cycles before changing its frequency by a factor of 2. This 
points out to a possible increase of the estimated S/N by a factor s/ri. The importance of accurate modelling of the 
expected waveforms is well appreciated. Presently, this is an important challenge for theorists (see, for example, [25]). 
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FIG. 3. The status of sensitivity of one of the LIGO interferometers 

Fig. 4 (taken from [18]) shows the accurately calculated S/N ratios for the initial laser interferometers. As expected, 
the S/N grows with the total mass M of the coalescing binary, because the signal gets larger. However, binaries with 
the total mass greater than about 80M Q radiate at too low frequencies before merging. At these low frequencies the 
detector noise is high, as seen in Fig. 3. Therefore, the S/N decreases again. Certainly, the equal mass binaries with 
the total mass greater than M ~ I0M Q can only be pairs of black holes, according to our present understanding of 
stars and stellar evolution. 

The LIGO-GEO collaboration has reported [26] the first upper limits on coalescing binary systems, as well as other 
possible g.w. sources. These limits were derived from observations at the currently achieved sensitivity. They are not 
yet significant from the astrophysical point of view, but, outside the frequency interval probed by the bar detectors, 
they are tighter than other upper limits experimentally established so far. The necessary steps for increasing the 
sensitivity in the initial and advanced ground-based interferometers are well recognised and are being taken up. 

The space-based interferometer LISA will be sensitive to g.w. in the interval I0~ 4 — 10°Hz, that is, to lower g.w. 
frequencies as compared with ground-based instruments. Some new types of astronomical sources will be accessible 
to this detector. LISA will consist of three spacecraft, forming an equilateral triangle of side 5 million km, in a 
heliocentric orbit, lagging behind the Earth by 20°. The phase shifts of the laser light traveing along all three sides 
of the triangle will monitor the light-travel distance between the small drag-free test masses inside the spacecraft. 
The nominal lifetime of the mission is 5 years. Fig. 5 shows the LISA design sensitivity curve and some of interesting 
sources of gravitational waves. This single figure attempts to show the noise curve of the entire assembly of spacecraft 
together with the strengths of g.w. sources of different nature - from quasi-periodic to stochastic - and therefore the 
figure should be treated with some care. The instrumental noise level is shown here in bins of 3 x 10~ s Hz, which are 
appropriate for a I year integration time. In other words, the r.m.s. instrumental noise */ S n (f) in units of (Hz)^ 1 / 2 
(see Eq.(14)) is multiplied with y3 x lQ~ s Hz at each frequency, thus producing a dimensionless spectral quantity. 
(For the latest amendments of the LISA noise curve see, for example, [27].) This is done mainly in order to emphasize 
that during this observational time the g.w. signals from many of the galactic white-dwarf binaries can be resolved 
and removed from the data. This means that their g.w. noise will not prevent us from seeing something much more 
interesting - a stochastic background of relic gravitational waves. The sharp drop of signal from the galactic binary 
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white dwarfs, shown in Fig. 5 at / = 2 x 10~ 3 Hz, illustrates this assumed operation with the data [18], [28], [29], and 
not the total lack of galactic binaries radiating at frequencies higher than 2 x 10~ 3 Hz. The continuation of this curve 
at / > 2 x 10~ 3 Hz shows the much smaller g.w. noise from unresolved extragalactic white dwarf binaries. In accord 
with this way of describing the dimensionless instrumental noise, the dimensionless g.w. signal amplitudes are also 
calculated in bins of 3 x 10~ 8 Hz around any given frequency /. In the next section, we will present more details on 
astrophysical sources for ground-based and space-based instruments. 
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FIG. 4. Signal to noise ratio in initial interferometers as a function of total mass for inspiral signal from binaries of equal 
masses and averaged over source inclination 



IV. GRAVITATIONAL WAVES AND ASTROPHYSICS 

It is common to divide the sources in groups, depending on whether they are accessible to ground-based or space- 
based instruments. It is also common to call them, respectively, the high-frequency and low-frequency sources. 



A. High-frequency sources 

We will focus on the sources that are likely to be the first sources detected by the ground-based instruments. 
It is clear from Eq. (10) that a powerful source of gravitational waves should involve large masses and relativistic 
velocities. A pair of neutron stars inspiraling toward each other under the influence of the g.w. radiation reaction 
force is a primary example. In the last few seconds of inspiral, before the neutron stars "touch" each other and 
merge, the frequency of gravitational radiation increases from 200Hz to about 1200Hz. The orbital velocity increases 
from v/c » 0.2 to v/c ~ 0.4. In the last few seconds the binary emits 1.1 x 10 53 ergs of energy in the form of 
gravitational waves. This is a 2% of its total rest-mass energy Mc 2 . What can be more powerful and efficient a 
source of gravitational radiation than this one ? Only a pair of even more massive and even more compact stars. 
Astronomers call them stellar-mass black holes. 

If the coalescence of a binary neutron star (NS+NS) or a binary black hole (BH+BH) happened in our Galaxy, it 
would be easily detectable by the already operating instruments. The problem is that these catastrophic events are 
expected to take place only once per a very long while. In order to have a reasonable chance of seeing, say, 3 events 
per year, we have to survey a large volume of space, which includes many galaxies. This means that the instrument's 
sensitivity should be so high, that the events could be seen from the edges of this large volume. This consideration 
explains the importance of theoretical evaluations of the event rates in a typical galaxy. 

The event rate of coalescing NS+NS systems is partially constrained by pulsar observations in our Galaxy. Three 
NS+NS binaries are known, each involving a pulsar, whose coalescence time is less than the Hubble time and is, on 
average, 3 x 10 8 years. Starting with these 3 binaries, one would evaluate the NS+NS event rate as 1 per 100 million 
years. However, we observe only about 1% of the galactic volume, so the coalescence rate can easily be raised to 
10 _6 ?/r _1 [30]. Very likely, the event rate for NS+NS systems is significantly higher than this estimate, if only because 
of the fact that not all NS+NS systems include a currently observable pulsar. The observational situation with black 
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holes is more uncertain. There are a dozen of BH candidates in X-ray binary systems, but they all arc in pairs with 
non-degenerate companions. So far, there is no observational evidence for NS+BH or BH+BH binaries. Nevertheless, 
one can make some evaluations on the basis of the star formation rates. 

It is believed that the neutron star progenitors have masses greater than 10M©, whereas the black hole progenitors 
have masses greater than 8OM . The Salpeter function for the star formation rate is 



dN / M N ; ! 



dtd(M/M Q ) \M Q 



yr 



Integrating this function over M and using the lower limit of integration, one finds the ratio of the expected numbers 
of progenitors: 



N(M > 80M Q ) _ /80M Q 
N(M > 10M Q ) ~ V1OM 



-1.35 

~ 0.06. 



It is reasonable to think that, despite all the complexities and differences in binary evolution, the ratio of coalescence 
rates will also be given by approximately the same quantity, 



K BH _ /8OM X " :r ' 







0.06. 



This expectation turns out to be in rough agreement with the results of detailed numerical population synthesis 
calculations. 

Numerical calculations take into account all the available observational information. Their advantage is in that one 
can follow not only the channels leading to the NS+NS, NS+BH, BH+BH systems, most interesting for gravitational- 
wave astronomy, but also other evolutionary outcomes, which allow comparison with observations in their own right. 
The population synthesis results cannot be less reliable than purely "observational" estimates, as they are controlled 
by the same available observational data. Discrepancies in final results exist because of uncertainties in astrophysics, 
not because one of the methods is inherently less reliable than another. This is especially true with regard to the 
NS+BH and BH+BH binaries, where the current purely "observational" evaluations would have to begin with zero. 

The results of conservative population synthesis calculations [31], [18] have been checked on their consistency with 
other evolutionary outcomes. The calculations show that the NS+NS rate is expected to be at the level 1Zns = 
3 x 10 _5 yr _1 , whereas the BH+BH rate is at least one order of magnitude lower. For further estimates we will take it 
at the level IZbh = 0.06 1Zns = 2 x 10 _6 j/r _1 . These rates for a typical galaxy, TZq, determine the rates for a given 
cosmological volume, IZy, which includes many galaxies. When deriving the 72.y, it is convenient to use a conservative 
estimate for the baryon content of the Universe. This brings us to the relationship 

n v w o.m G 



1 Mpc 



Thus, it is expected that within the volume of radius r = 100 Mpc and during 1 year, there will be 3 of NS+NS events 
and only 0.2 of NS+BH or BH+BH events. The increase of the radius to r = 200 Mpc increases the volume and the 
event rates by a factor 8. It is interesting to note that, despite all the diversity of approaches in the literature, there 
exists some tendency to convergence of the final results [32] , [33] , [34] , [35] , [36] . 

The derived event rates are the basis for the calculation of the expected detection rates. The important fact is 
that the mass of a typical neutron star is 1.4 M Q , whereas the mass of a typical black hole is (10 — 15) M Q . The 
averaged mass of the observed black hole candidates in the X-ray binaries is Mbh — 8.5M Q . A pair of black holes is 
a more powerful source of gravitational waves than a pair of neutron stars, and therefore black hole binaries can be 
seen by a given instrument from much larger distances. When it comes to the detection rates, the lower event rate of 
BH+BH sources is more than compensated by their larger masses. Indeed, the optimal signal to noise ratio is [15], 
[37], [18]: S/N cx M 5 ^ 6 /r. At a fixed S/N, the detection volume is proportional to r 3 and therefore to M 5 ^ 2 . The 
detection rate V for binaries of a given class is the product of their event rate IZy and the detector's registration 
volume cx M 5 ^ 2 for these binaries. Therefore, one obtains 

Vbb = Ubh (Mbh \ 5/2 = QM (S_5M^^ ^ 
V NS 1l N s \Mns J \IAMq) 

This remarkable result is consistent with the more accurate numerical calculation of S/N displayed in Fig. 4. 
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It is seen from Fig. 4 that g.w. signals from NS+NS at r = lOOMpc cannot be regarded detectable by initial 
interferometers. The situation is much better for havier pairs of NS+BH and BH+BH. If the total mass of a BH+BH 
binary is near (20 — 30)M & , then S/N w 2 even if the binary is placed at r = 200Mpc. It follows from the event rate 
IZbh discussed above, that in this larger volume one expects a couple of BH+BH events per year. The simultaneous 
observations on two or three instruments will significantly diminish the probability of false alarms to such events. 
This is why it is argued in Ref. [18] that the coalescing black holes will probably be the first sources detected by the 
initial ground-based interferometers, when they reach their planned sensitivity. Of course, the discussed estimates 
are statistical by the very nature of things, and they have significant systematic unceratinties. It will not be very 
surprising if the reality is somewhat better or somewhat worse than what the mean values suggest. It is important, 
however, that even the most pessimistic evaluations of the NS+NS and BH+BH rates indicate that there should be 
many detections per year in the advanced interferometers. 

If the coalescing black holes are detected first, it is likely that at the beginning we will only have a proof that the 
objects are black holes in astronomical sense - an inspiraling pair of heavy compact masses. The next step will be 
to try to understand their real physical nature. The general-relativistic black holes possess event horizons, which are 
supposed to merge into the resulting black hole, which will then emit a damped train of ringdown waves at specific 
frequencies, and so on [37], [38]. This fascinating physics will be testable when the good quality data are available. 
This analysis will also require the continuation of the intense effort on the side of analytical and numerical calculations 
(for a recent review of numerical relativity, see [39] ) . 

Direct detection of first sources will be not the end, but only the beginning of the observational g.w. astronomy. 
In the long run, the aim of g.w. science is to explore a great variety of sources, many of which can hardly be seen in 
electromagnetic radiation. Needless to say that there is also a great chance of discovering new and totally unexpected 
sources. In addition to coalescing binary stars, many other sources will eventually be detected. For instance, the 
long-recognised importance of the core collapse of massive stars [40], [41] as g.w. sources has been reinforced by the 
mounting evidence of asymmetries during the supernova explosions, by the likelihood of forming and quick collapse 
of very massive early stars, by the association of supernovae events with gamma-ray bursts, etc. (For a recent review 
and extensive list of references, see for example [42].) The tidal disruption of a neutron star by its companion, the 
various sorts of stellar instabilities [43], the slightly deformed spinning neutron stars, young pulsars, and low-mass 
X-ray binaries [44], [45] - are also the astrophysically important and interesting g.w. sources that will be studied, very 
likely, by advanced detectors. 



It is common to call the low-frequency g.w. sources as "sources for LISA". Some of them are displayed in Fig. 5. 
As explained in the previous section, the dashed line shows the g.w. confusion noise from the binary white dwarfs, 
WD+WD, mostly concentrated in the disk of our Galaxy. LISA will not only detect thousands of WD+WD systems 
radiating at / > 2 x 10~ 3 Hz, but is capable of doing this so accurately that their contributions can be removed 
individualy from the data. Some known binaries consisting of degenerate and normal stars will also be detectable. 
Their angular coordinates and distance can be measured with high precision [46]. In fact, the well identified galactic 
binaries serve as guaranteed sources for LISA, and they will help to test LISA's performance. 



B. Low- frequency sources 



le-18 




FIG. 5. Some interesting sources in comparison with LISA sensitivity 
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There is growing evidence for the existence of binary supermassive black holes (SMBH, M > 1O 6 M0) in the centers 
of merging galaxies. (For a recent review of SMBH formation see [47]). Certainly, a coalescing pair of SMBH is 
an extremely bright g.w. emitter. These sources are located at cosmological distances, so that the red shift of the 
incoming gravitational waves becomes important. But if the total mass of the pair is not significantly larger than 
10 7 M Q , the coalescing pair can still be visible by LISA at frequencies up to / w 10~ 4 Hz, even if the source is located 
at the red shift z = 3. The upper line in Fig. 5 shows the effective gravitational wave amplitude, relative to the plotted 
LISA instrumental noise, for a source consisting of two SMBH. The dot shows the signal when the inspiraling pair 
was only 1 year prior to merger. The presence of spins of SMBH changes the LISA signal [48]. The trouble with 
these super-powerful sources is that their event rate is very uncertain. One of difficult issues is whether a SMBH 
forms by direct collapse of gas in deep galactic potential wells or by hierarchical build-up of pre-galactic structures. 
Nevertheless, it is estimated that LISA may see 0.1-1 events per year, or maybe a factor of 10 more [49]. 

It seems certain that, from time to time, a SMBH will be closely approached by nearby compact stars - white 
dwarfs, neutron stars and stellar-mass black holes. A compact star orbiting a SMBH is a powerful g.w. source. 
However its orbit and the emitted waveform depend strongly on the spin of SMBH. Fig. 5 shows two curves (taken 
from [19]) describing the performance of a 1OM BH completing its inspiral toward a 10 6 M Q SMBH. The left curve 
corresponds to a non-spinning SMBH and a circular orbit of a stellar-mass BH. The right curve describes the signal 
from a prograde, circular, equatorial orbit around a nearly maximally spinning Kerr SMBH. Both curves begin at 
frequencies when the time to merger is 1 year. A serious theoretical problem is the construction of reliable templates 
for these sources. One should expect that a typical orbit will be eccentric, non-equatorial, and subject to the radiation- 
reaction force corrections. Correspondingly, the waveforms could be extremely complicated. Some progress in this 
area is reported in [50]. 

Finally, Fig. 5 shows the expected level of relic gravitational waves. This is a fundamentally important signal from 
the very early Universe. Its explanation requires cosmological notions and some elements of quantum physics. Since 
the relic gravitational waves and primordial density perturbations is presently one of the most fascinating and active 
areas of research, we devote to it a separate section. 



V. GRAVITATIONAL WAVES AND COSMOLOGY 



A. Generation of relic gravitational waves and primordial density perturbations 

In many situations one can neglect the non-linearity of the gravitational field, that is, the interaction of gravitational 
waves with other gravitational fields and with themselves. However, this is not always the case. The most dramatic 
example is the interaction of gravitational waves with the strong variable gravitational field of the very early Universe. 
A gravitational wave can be thought of as a harmonic oscillator, while the smooth variable gravitational field of the 
surrounding Universe as a gravitational pump field. The g.w. oscillator is parametrically coupled to the gravitational 
pump field. This specific coupling is a consequence of the non-linear nature of the Einstein equations. The coupling 
provides a mechanism for the superadiabatic (parametric) amplification of classical waves and the quantum-mechanical 
generation of waves from their zero-point quantum oscillations [51]. The word "superadiabatic" emphasizes the fact 
that this effect takes place over and above whetever happens to the wave during very slow (adiabatic) changes of 
the pump field. That is, we are interested in the increase of occupation numbers, rather than in the gradual shift of 
energy levels. The word "parametric" emphasizes the mathematical structure of the wave equation. It is a change of 
a parameter of the oscillator caused by the pump field, namely, a sufficiently rapid variation of its frequency, that is 
responsible for the considerable increase of energy of that oscillator. 

It is common to write the perturbed gravitational field of a homogeneous isotropic universe in the form: 



ds 2 = a 2 (77)[-d77 2 + (Sij + h^dx'dx^}. 



(15) 



The gravitational field perturbations hij(r],x) can be expanded over spatial Fourier harmonics e , where n is a 
constant (time-independent) wave-vector, 



h n (v)e lr 



(16) 



The polarisation tensors Pjj(n), s = 1, 2 have different forms, depending on whether the hij represent gravitational 

s 

waves or density perturbations. In the case of gravitational waves, the Pij describe the two familiar "plus" and "cross" 
polarisations introduced in Sec. 2. In the case of density perturbations, the polarisation tensors are: 
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i 12 2 /-riin-i 1 

Pij(n) = yj -Sij, P ij (n) = -V3-^ + -j=5 ij . (17) 

The Einstein equations for the gravitational field perturbations with the polarisation tensors (17) can only be 
satisfied if the hij are accompanied by perturbations in the density of matter. This is why this class of perturbations 

s 

is called density perturbations. The difference between Pij(n) for, respectively, gravitational waves and density 
perturbations is responsible for the difference in polarisation patterns of the CMB radiation, caused by these two 
classes of gravitational perturbations (see Sec. 5.2 below). 

s st 

For a classical field hij, the quantities c n , c n are complex numbers. For a quantized field, they arc annihilation and 
creation operators satisfying the conditions 

f c n, c m ] = S s , s S 3 (n - m) , c n |0) = , 

where |0) (for each mode n and s) is the initial vacuum state defined at some 770 in the very distant past, long 
before the superadiabatic regime for the given mode has started. The normalization constant C is determined by the 
requirement that initially each mode contained only the zero-point energy \hio. Then, C = y/Wnlpi for gravitational 
waves and C = V^i-Klpi for density perturbations, where Ipi — (Gh/c 3 ) 1 / 2 is the Planck length. Obviously, the initial 
vacuum amplitude, and the entire field, should vanish, if the Planck constant U is formally sent to zero. 

The calculation of quantum-mechanical expectation values and correlation functions provides the link between 
quantum mechanics and macroscopic physics. Using the representation (16) and definitions above, one finds the 
variance of the gravitational field perturbations: 

(0|Mf,,x)fc«(„,x)|0) = 0-2 / ^ E IMtfl 2 — ■ (18) 

Z7T Jo s =l,2 71 

The quantity 

h 2 (n, V ) = ^n 2 J2 \Uv)\ 2 (19) 

8=1,2 

gives the mean-square value of the metric (gravitational field) perturbations in a logarithmic interval of n and is called 
the (dimensionless) power spectrum. The power spectrum of metric perturbations is a quantity of great observational 
importance. It defines the temporal structure and amplitudes of the g.w. signal in the frequency bands of direct 
experimental searches. It is also crucial for calculations of anisotropy and polarisation induced in CMB by relic 
gravitational waves and by other gravitational field perturations. 

To find the power spectrum at any given moment of time (for instance, today or at the moment of decoupling of 

s 

CMB from the rest of matter) we need to know the mode functions h n (v) at those moments of time. In the case 

s s s 

of gravitational waves, the mode functions M n (»?) (where M n (ry) = a(ri)hn(v)) are governed by the equation for the 
parametrically disturbed oscillator [51]: 



s" s 



n 

a 



= . (20) 



The equation describing density perturbations in the very early Universe can also be reduced to the equation very 
similar to Eq. (20), and the generating mechanism will work without change. As soon as the pump field (represented 
by the cosmological scale factor a(r))) is known, and since the initial conditions arc fully determined, the mode functions 

8 

hn(v)i as wen as other properties of the generated fields, are unambiguously calculable. The important (and, strictly 
speaking, unknown) era of cosmological evolution is the stage preceeding the radiation-dominated era. We call it an 
initial (i) era and characterize, for simplicity of calculations, by a set of power-law scale factors: a(rj) — l \ri\ 1+f3 , 
where l Q and f3 are constants. 

The main properties of cosmological perturbations generated as a result of superadiabatic (parametric) amplification 
of their zero-point quantum oscillations are as follows (for more details see [18], [52] and references there): 

1. The initial vacuum state is described by a Gaussian wavefunction. As a result of quantum- mechanical Schrodinger 
evolution, the vacuum state transforms into a multi-particle state known as a squeezed vacuum state. The distributions 
of amplitudes and phases acquire strongly unequal variances. In the general expression for the gravitational field mode, 
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h n = A\ sin(ra? + fa) cosn • x + A 2 sm(nr/ + fa) sinn • x, 



the amplitudes A\ and A 2 are drawn from a broad Gaussian distribution, whereas the phases fa and fa are practically 
fixed and equal up to ±ir. The field is a stochastic collection of standing waves and is characterised by a strongly 
modulated power spectrum. This complicated statistical picture of generated cosmological perturbations is often 
replaced in the literature by a single word: "Gaussian" . 

2. The generated gravitational field perturbations act on all sorts of matter together. There is no reason why the 
inhomogeneitics in different sorts of matter (if more than one component of matter was dynamically important in the 
very early Universe) should be displaced and move with respect to each other, which would constitute the so-called 
isocurvature, or entropy, perturbations. This is why the generated density perturbations are often called "adiabatic". 

3. The primordial spectrum (i.e. the spectrum before processing at the radiation-dominated and matter-dominated 
stages) as a function of the wave-number n is fully determined by the variable pump field as a function of time rj. Every 
interval of spectrum that can be meaningfully approximated as a power-law function of n, was generated by an interval 
of a power-law evolution a(r]) cx |ry| 1+ ' 3 . Concretely, for the spectrum of Eq. (19) one finds: h 2 (n) cx (lpi/l ) 2 n 2 ^ +2 \ 
Specifically for density perturbations, one often uses the spectral index n, related to (5 by n = 2/3 + 5. If the very early 
Universe was goverened by a scalar field (the central assumption of inflationary scenaria), then, at every power-law 
interval of evolution and, hence, at every power-law interval of the generated spectrum, there must be (3 < —2 and 
n < 1. The "red" spectra (3 < —2 (n < 1) possess a serious theoretical difficulty: the mean-square value of the field 
becomes power-law divergent in the limit of very long waves (lower limit of integration in Eq.(18)). The case (3 = — 2 
(n = 1) is called the flat, or Harrison- Zeldovich-Peebles, or "scale-invariant", spectrum. 

4. Gravitational waves are generated inevitably, wheras the generation of density perturbations requires additional 
assumptions about the coupling of matter fields to gravity. The primordial (unprocessed) amplitudes of density 
perturbations can be as large as g.w. amplitudes, but never much larger. The observed CMB anisotropy in lower 
multipoles (caused by primordial gravitational field perturbations) may have nothing to do with quantum mechanics, 
but if it does, the contribution of relic gravitational waves to lower multipoles must be substantial. 

5. The parametric mechanism is universal, and the generation of primordial cosmological perturbations takes place 
regardless of whether the Universe 13 billion years later will appear to astronomers as spatially- flat, or not. In any 
case, if Qtotai is not identically 1, the present-day Qtotai is usually regulated "by hand" through the duration of the 
initial era. This parameter of duration does not affect seriously the amplitudes and spectral slopes of primordial 
perturbations. 



The direct detection of relic gravitational waves will be fundamentally important for uncovering the physics of 
the very early Universe. In Fig. 6 we show the expected spectrum of today's r.m.s. amplitudes h(v) (square root 
of Eq. (19)) as a function of frequency v in Hz. The graph shows the piece-wise envelope of the spectrum and 
ignores its oscillations. Almost everything in this graph is the processed spectrum; the primordial part survives only 
at frequencies v < i/h, where vh = c/Xh = H w 2 x 10~ 18 Hz is the Hubble frequency. This particular spectrum 
was derived under the assumption that a significant fraction of the observed large-scale CMB anisotropy is caused by 
relic gravitational waves and that the primordial spectral index is = — 1.9 (n = 1.2). This value of n follows from 
the COBE data [53] and it has been recently reinforced, albeit with broad error bars, by more sophisticated analysis 
[54]. The evaluations of n based on larger data sets usually lead to smaller n's, varying in a narrow interval around 
n = 1. However, this precision appears to be artificial (caused by the excessively rigid apriori assumptions about the 
tested models) and the story does not seem to be over. In any case, we use Fig. 6 for the analysis of detectability of 
relic gravitational waves. 

The r.m.s. values h{v) are directly entering the detectability evaluation, but they also determine the Q. gw {v)- 
parameter, which is useful for comparison of the g.w. background with other energy components. This parameter 
can be calculated according to the formula 



For example, one has h{u) w lO" 205 , Q.{v) w lO" 11 at v = 10~ 3 Hz, and h{v) w 10" 25 , Q{v) w 1CT 10 at v = 10 2 Hz. 
It is seen from Fig. 5 that the part of the spectrum accessible to LISA is higher than the instrumental noise and can 
be measured. The ground-based advanced interferometers are also promising. The target value of the LIGO-II at 
v = 10 2 Hz is h ex — 1(U 23 . The gap in two orders of magnitude can be covered by cross-correlation of the outputs of 



B. Detection of relic gravitational waves 
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two or more detectors. The S/N will be better than 1 if the common integration time exceeds 10 6 sec. This does not 
seem to be a hopeless task. 

It is possible that relic gravitational waves will be first observed indirectly, with the help of polarisation measure- 
ments of CMB. The problem is to distinguish the polarisation pattern caused by gravitational waves from that caused 
by density perturbations. The polarisation arises as a result of the Thompson scattering of CMB photons on free elec- 
trons [55]. To produce a net polarisation, the electrons should be illuminated by CMBR having a non-zero quadrupolc 
anisotropy. The two polarisation patterns are distinguishable if the quadrupole anisotropics are distinguishable. For 
our purposes of illustration, it is sufficient to consider the emission (rather than scattering) of electromagnetic waves 
by free electrons. The electrons are set in motion by, respectively, gravitational waves and density perturbations. 



log h(v) 




logv, Hz 

FIG. 6. Expected envelope of the spectrum h(y) for the case f3 = —1.9 (n = 1.2) 

The motion of free particles in the field of a linearly-polarised gravitational wave is shown in Fig.l. Imagine that the 
moving particles are free electrons in the early Universe, rather than free mirrors of an interferometer in laboratory, 
that we were discussing in Sec. 3. Then, the directions of the induced osillations of the electrons, indicated by arrows 
in Fig.l, are, at the same time, the directions of the electric fields of the electromagnetic waves emitted by these 
oscillating electrons. The pattern of arrows seen on this figure is the pattern of polarisation components that will be 
seen on the sky. 
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FIG. 7. Motion of free particles in the gravitational field of a density perturbation 

The motion of free electrons induced by the gravitational field of a density perturbation is different. Two polarisation 
tensors of a density perturbation are given by Eq. (17). The first one describes a trivial purely spherically-symmetric 
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deformation of a sphere of free particles. The second one, subject of our interest, describes the deformation which has 
the angular dependence of the spherical harmonic Yi m (9, </>) with 1 = 2 and m = 0, in contrast to the combination of 
Yi m (0, <fr) with I = 2 and m = ±2, attributable to a gravitational wave. In both cases, the polar axis z is taken along 
the wave-vector n of the perturbation. Fig. 7 shows a ring of electrons set in motion by the gravitational field of a 
density mode. The pattern of these arrows is, at the same time, the pattern of polarisation components generated by 
a density perturbation. Comparing Fig.l and Fig. 7 one concludes that the polarisation patterns are distinguishable, 
if one can study the distrubution of polarisation Stokes parameters over a sufficiently large portion of the sky. In 
real conditions, the free electrons will be influenced by the superposition of many perturbation modes with arbitrary 
wave-vectors, but the net difference between these two sorts of perturbations should survive. 

The gap between the expected g.w. signal and the polarisation detection capabilities is relatively small. Taking 
into account the current impressive activity in this area, one can hope that some decisive observational information 
about relic gravitational waves can be obtained pretty soon. 



C. Gravitational waves and inflation 



The theoretical and experimental studies of relic gravitational waves are endangered by absurd claims that are 
prevalent in inflationary literature. Inflationists claim that the amount of relic gravitational waves should be zero if 
the primordial spectrum is flat, that is, if n = 1. There seems to be no need to bother about relic gravitational waves, 
as the observations indicate that n may indeed be close to 1. Since the "pillars of inflation" are popular in a part of 
astrophysics community, and sometimes are said to be "confirmed", it is important to put matters straight. 

Inflationary scenario operates with a scalar field p> and the scalar field potential V(ip). Having accepted the 
general concept of parametrically amplified quantum fluctuations, inflationary theorists are performing their own 
calculations. The quantum-mechanical content of these calculations is usually limited to vague words, such as that 
"inflation amplified quantum fluctuations onto macroscopic scales" . Being unsure why and where the Planck constant 
H should enter the calculations, inflationists never write it explicitcly; and when it is written implicitely, in the form 
of the Planck mass, Mpi = (Tic/G) 1 / 2 , it always stands in the wrong place, in the denominator of the final expression 
instead of the nominator. With this sort of "quantisation" , inflationary theorists derived their contribution to the 
subject of cosmological perturbations - the "standard inflationary result" . The "standard inflationary result" predicts 
the infinitely large amplitudes of today's density perturbations in the limit of the flat spectrum n = 1 . Indeed, one will 
always be able to recognise in inflationary papers the evaluations relating the final (/) amplitudes of the perturbations 
to the initial (i) values of (p and other quantities: (Sp/p)f ~ (hs)f ~ (Of ~ (C)i ~ (H 2 /(p)i ~ (V 3/>2 (p)/V' (<p))% ~ 
Hi/y/1 — n. The denominator of the last expression is zero for n = 1. The nominator is the Hubble paparameter 
Hi at the initial stage. Hi is much larger than the Hubble parameter at the subsequent radiation-dominated stage, 
and, in any case, Hi is not zero. Therefore, the predicted final amplitudes go to infinity, if the spectral index n = 1. 
Inflationists are hiding this absurd prediction of infinitely large density perturbations by composing the ratio of the 
gravitational wave amplitude hx to the predicted divergent amplitude of the scalar metric perturbations hs (the so 
called "tensor-to-scalar ratio" or "consistency relation": hx/hs ~ \/l — n) and declaring that it is the amount of 
gravitational waves that should be zero, or almost zero, at cosmological scales and, hence, down to laboratory scales. 

Certainly, the "standard inflationary result" is in full disagreement not only with the theoretical quantum mechanics, 
but with available observations too, as long as the error-boxes of the observationally derived spectral index n are 
centered at n w 1 and include n = 1. To be consistent with inflationary predictions, the data should not allow "blue" 
spectra n > 1, as the scalar field cannot produce them, and the density amplitudes should go to infinity, when one 
processes the data assuming that n = 1. This spectacular failure of inflationary calculations is systematically painted 
by inflationists and their followers as a great success. The most recent example is the analysis of WMAP data [56]. 
The authors praise and follow inflationary derivations, and conclude that the "tensor/scalar ratio r" is "consistent 
with zero". The "standard inflationary result" is written in that paper (their formula (17)) as: 

2 _ V/MU 
Ak _ 24^ ey ' 

and the "tensor/scalar ratio r" (their formula (18)) as: 

r = 16ey. 

Combining the second formula with the first one, one can easily see that if the WMAP data demonstrate that "r is 
consistent with zero" , then the WMAP data should also be consistent with an infinite numerical value of the density 
amplitudes A^j and, hence, with an infinite numerical value of the induced CMB anisotropies. If the WMAP data are 
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not consistent with such an infinite numerical value of density perturbations, then the only "implication for inflation" 
that follows from the WMAP observations is that the single concrete formula derived by inflationists - their "standard 
inflationary result" - is shown to be wrong. 

It seems to the author that overenthusiasm for inflation has reached unscientific, even ecclesiastical proportions. 
For example, this is how inflation is characterized in the educational programm "Astronomy" of the Smithsonian 
Institutions [57]: "...the inflationary scenario is the best current theory of the Universe... It has met four critical obser- 
vational tests..." . It is unclear which 4 tests the inflationary school credits to itself, rather than to direct consequences 
of quantum mechanics and general relativity, but one can recall that even general relativity was characterized until 
quite recently as a theory that has met only 3 critical observational tests (1 of which, gravitational red-shift, is not 
a test of specifically general relativity). As for professional papers, one reads almost every day claims that the CMB 
and galaxy surveys are "in spectacular agreement with an inflationary A-dominated cold dark matter cosmology" 
(compare, for example, with [58] and [59]). Somewhere in the text, authors usually admit that, say, the observed 
quadrupole anisotropy is way out of the predicted value, and that the probability of finding such a result within the 
"standard" model is 1.5 x 10~ 3 [60], [58]. Surely a theory which is admitted to have only a one on a thousand chance 
of being consistent with one of its crucial observational tests is not in "spectacular agreement" with the cosmos we 
are trying to understand and should not be a subject given to self-congratulation. 

D. Gravitational waves and quadrupole anisotropy 

The accurate measurements of CMB by WMAP reiterate the issue of the gravitational wave contribution to the 
lower order multipoles. The best strategy is to rely on conclusions of general physics and to use the minimum number 
of extra hypotheses. If the general considerations suggest (see Sec. 5.1) that the contributions of gravitational waves 
and density perturbations should be of the same order of magnitude, it is this conclusion that should be tested most 
thoroughly. 




FIG. 8. The WMAP data and some theoretical models 



Fig. 8 shows [61] the WMAP observational points (triangles) and the best fit curve (solid line) of the A-dominated 
cosmology without gravitational waves. The predicted quadrupole (I = 2) anisotropy is a factor of 8 higher than 
the actually observed value. The thin solid line (reproduced from [52], fig. 12) shows the contribution of density 
perturbations alone in a model with Zdec — 

1000, z eq = 5000, n = 1 that fits the position and value of the peak at I = 
220. This line is surrounded by the la uncertainty belt (shown by two dashed lines) arising due to the lack of ergodicity 
on a 2-sphere. At sufficiently large I's the belt is approximately symmetric and its size is AC; w y/2/(2l + l)Cj. But 
the asymmetry grows towards small i's, and specifically at I = 2 the size of the deviation down is only 0.4 part of the 
deviation up [62] . [The accurate displaying of this fundamental uncertainty makes the actual quadrupole even further 
out of the uncertainty belt of the A-dominated cosmology, than what is implied by the usually plotted symmetric 
"cosmic variance" .] Comparing the thin line with the data points, it is difficult to avoid the conclusion (advocated 
in [52]) that in fact there exists an excess, rather than a deficit, of power at small multipoles, and the most natural 
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explanation of this excess is the anticipated contribution of gravitational waves. The increase of the spectral index 
up to n = 1.2 makes the agreement with the observed quadrupole even better and implies a somewhat larger amount 
of gravitational waves [52]. One should remember, however, that all the experimental C; data, together with all the 
cosmological parameters, is a finite set of numbers. At the same time, in our hands is the (strictly speaking, unknown) 
shape of the primordial spectrum, i.e. a continuous function and an infinite set of numbers. A perfect agreement with 
any observed C;'s and practically any set of cosmological parameters, not only A = 0, can be achieved at the expense 
of the properly chosen primordial spectrum. Unfortunately, the era of "precision cosmology" is still at some distance 
from us. 

The nature of the observed quadrupole anisotropy deserves special attention. Most likely, it is caused by superpo- 
sition of very long gravitational and density waves. It is known [63] that a gravitational wave produces the Y~ 2 , 2 and 
y 2 ,-2 CMB anisotropy, whereas a density perturbation produces the F 2 .o CMB anisotropy. We have discussed this 
distinction in Sec. 5. 2 in connection with the CMB polarisation. The actual quadrupole distribution over the sky was 
measured by COBE [64]. In Galactic coordinates, 

Q{0,<t>) = Qi(3cos 2 - l)/2 + Q 2 sin20cos0 + Q 3 sin 20 sine/) + Q 4 sin 2 cos 2(f) + Q 5 sin 2 sin 20, 

where the measured (least noisy) components are [65]: 

Qi = 19.0 ±7.4, Q 2 =2.1±2.5, Q 3 = 8.9±2.0, Q 4 = -10.4 ± 8.0, Q 5 = 11.7 ±7.3. (21) 

Even if this Q(9, 4>) is produced by a single gravitational wave or a single density perturbation, it is only in a special 
coordinate system that it can be reduced to the combination of Y" 2i2 and ¥2,-2 or to Y" 2 .o, respectively. To find out 
what we are dealing with, we have to build invariants, that is, quantities independent of the rotation of the coordinate 
system. One of invariants is 

Q 2 rms = (4/15)[(3/4)Q 2 + Q\ + Ql + Ql + Q\], 

another one (see, for example, [66]) is 

D = (4/5 3 / 2 )[(l/4)g!(Q 2 + Q 2 /2 + Q 2 /2 - Q\- Q 2 ) + 2Q 2 Q 3 Q 5 + Q 4 (Q 2 2 - Q§)]. 

Qrms is always positive, wheras D can be negative, but it always satisfies the condition \D\ < Qrms- For a pure 
density perturbation, \D\ = Qf ms ; and f° r a pure gravitational wave, D = 0. 
Calculating the invariants and the formal errors from the data set (21) we find 

Q rms = (12.6 ± 3A)fj,K, - D 1/3 = (6.9 ± 12.9)fiK 

As expected, the available noisy data do not allow one to prefer one of hypotheses over another. But, for sure, 
there is no indications whatsoever that the quadrupole anisotropy is produced by a density perturbation alone. If 
anything, the central values of Q rm s and D indicate that the contribution of gravitational waves should be substantial. 
Hopefully, the WMAP quadrupole data will be more accurate, and then this analysis should be repeated. 

VI. SUMMARY 

Gravitational-wave physics is a mature and at the same time a very young science. In a sense, the relativistic gravity 
(general relativity) itself is still a young science. The enormous progress in technical developments and observational 
verifications of the theory is accompanied by difficult issues of its adequate description and interpretation, the necessity 
of bringing it to a closer contact with other branches of physics. 

It is rumored that the Nobel recognition eluded the great astronomer E. Hubble because of his reluctance to say 
about his discovery what the establishment wanted him to say. Apparently, the scientific integrity of E. Hubble 
allowed him to say what he believed he discovered - the nonstationarity of the system of nearby galaxies - whereas 
he was required to admit that he discovered the "expansion of space". But the "expansion of space" is still alive 
and well. [For example, "As bizarre as it may seem, space itself is expanding - specifically, the vast regions of space 
between galaxies" [57].] It is regularly proposed to be measured. The logic seems to be impecable. If "space" expands 
by a factor of 2 in 10 billion years, why would not the Earth or an atom expand by 10% in 1 billion years ? The 
gravitational-wave research is plagued in a similar fashion. It is often stated that "gravitational waves are oscillations 
of space-time itself . The next phrase seems to be logically unavoidable: "gravitational waves act tidally, stretching 
and squeezing any object that they pass through" . If this phrase were correct, we would never be able to notice 
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gravitational waves. The device measuring, say, the displacements of free mirrors in an interferometer would be 
"stretched and squeezed" as well. In this situation, we can probably find comfort in the wise observation [67]: "I 
agree that much of what one reads in the literature is absurd. Often it is a result of bad writing, rather than bad 
physics. I often find that people who say silly things actually do correct calculations, but are careless in what they 
say about them." 

It seems to the author that, in the long perspective, the value of the gravitational wave research will be in its 
influence on the fundamental physics. Meanwhile, let us hope that the next gravitational- wave update will be devoted 
to the fascinating nature of concrete astrophysical sources of gravitational waves detected by the existing instruments. 
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